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Abstract 

The status of dark matter searches with inorganic scintillator detectors at 
Boulby mine is reviewed and the results of tests with a CsI(Tl) crystal are pre- 
sented. The objectives of the latter experiment were to study anomalous fast 
events previously observed and to identify ways to remove this background. Clear 
indications were found that these events were due to surface contamination of 
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crystals by alphas, probably from radon decay. A new array of unencapsulated 
Nal(Tl) crystals immersed either in liquid paraffin or pure nitrogen atmosphere 
is under construction at Boulby. Such an approach allows complete control of the 
surface of the crystals and the ability to remove any surface contamination. First 
data from the unencapsulated Nal(Tl) do not show the presence of anomalous 
fast events. 
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1. Introduction 



The UK Dark Matter Collaboration (UKDMC) has been operating encapsulated 
Nal(Tl) detectors at the Boulby Mine underground site for several years [[[]. Com- 
petitive limits on the flux of weakly interacting massive particles (WIMPs), that may 
constitute up to 90% of mass in the Galaxy, have been set by one of these detectors 
using pulse shape analysis (PSA) to distinguish between scintillation arising from back- 
ground electron recoils and that due to nuclear recoils [0, |3|. Discrimination is possible 
because the sodium and iodine recoils expected from elastic scattering by WIMPs have 
faster mean pulse decay time than for electrons |Q. Traditionally, because Nal is hy- 
groscopic, detectors are fabricated using an outer copper encapsulation with glued-in 
quartz windows plus additional thick (typically > 100 mm) quartz light guides to shield 
the crystal from photomultiplier activity. However, this design limits detector sensitiv- 
ity because it prevents access to potential background sources on the crystal surfaces. 
The importance of Nal(Tl) surfaces has been highlighted recently by indications that 
they might be a source of anomalous fast time constant events seen so far at similar 
rates in many dark matter experiments with encapsulated Nal(Tl) crystals JT], [|, |^. [7J. 
Greater access would allow improved control of potential contaminants there and hence 
a possible reduction in such events, leading to greater sensitivity to WIMPs. 

In this paper we report the results of a study of anomalous fast events with Nal(Tl) 
and CsI(Tl) detectors and a recipe to suppress their rate. We confirm also that unen- 
capsulated Nal(Tl) crystals normally do not have this background. 

2. Anomalous fast events 

The observation of anomalous fast events in the UKDMC encapsulated Nal(Tl) 
detectors was first reported by Smith et al. [|TJ. These events are faster than typical 
electron recoil pulses and even faster than nuclear recoil pulses 0. Figure 0a shows 
typical time constant distribution of events from one run with 330 kg x days exposure 
of a 5.2 kg Nal(Tl) encapsulated crystal. The time constant distribution of events 
collected during calibration runs, when the crystal was irradiated by photons from a 
60 Co source, is shown in Figure |l]b for comparison. Both distributions are fitted to a 
log(Gauss) function, shown as a solid line in Figure [I]. The full data analysis procedure 
is described elsewhere 0, || [5J [UJ. A "bump" of anomalous fast events is clearly 
seen on the left edge of the distribution shown in Figure |I|a. Note the absence of the 
bump in Figure 0b, which implies that the bump is not due to any non-uniformity of 
the crystal. 

There is an excess of observed events over the log(Gauss) fit also at high values 
of time constant (see Figure p. This can be explained assuming stochastic pile- up of 
single thermal photoelectrons |TD[ , occasional afterpulses and fluorescence of the crystal 
after the scintillation pulse. The effect is seen for both "data" and "calibration" runs 
and does not interfere with the search for events faster than electron recoil events, such 
as nuclear recoils from WIMP interactions with matter or other kinds of fast events. 

Figure [2] show time constant distributions for another run with the same encapsu- 
lated crystal (a) together with the results of a calibration run with a neutron source 
(b). Both neutron- induced and gamma- induced events are seen on the distribution 
plotted on Figure 0b. The fits to gamma-induced events (electron recoils) are shown 
by dashed curves. The fits to anomalous fast events (a) and neutron-induced events 
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(nuclear recoils) (b) are shown by dotted curves. From the comparison of Figures 
0a and 0b we can conclude that anomalous fast events are faster than nuclear recoil 
events, expected from WIMP-nucleus interactions, and cannot be explained by WIMPs 
or neutron background. 

Similar fast events with comparable rate have been seen also by the Saclay group || . 
One of the Saclay crystals |TTJ, of size, growth technique, manufacturer, and housing 



similar to those used by the DAMA collaboration |L2|, [L3[, has been moved to Boulby 
(as a result of a collaboration between UKDMC and Saclay) and is currently collecting 
data. We confirm the presence of the population of fast events in this crystal with a 
rate similar to that seen in the UKDMC detectors (see also [|ll|, 14] for discussion and 
Figure |[) . 

Smith et al. suggested that the anomalous fast events could be due to MeV 
alphas. To account for the rate at low energies (10-100 keV) the alphas would need 
to deposit a small fraction of their initial energy at the crystal surface. Intrinsic bulk 
contamination of the crystal by uranium and thorium (measured to be at the level 
of about 0.1 ppb) is certainly not enough to explain the observed high rate at low 
energies. External incoming alphas from surrounding materials (PTFE - polytetraflu- 
oroethylene, quartz windows) cannot easily explain the observed spectrum: fine tuning 
of model parameters, such as a dead layer of scintillator, is needed and a very high 
contamination of the material by uranium or thorium (about 1 ppm) is required as 
well. Moreover, the time constant of the incoming alphas is not matched well to that 
of the fast events 0. 

Intrinsic surface contamination of the crystal by an alpha-emitting isotope has 
recently been discussed as a source of the anomalous fast events J7|. Recoiling nuclei 
from radon decay can be implanted into the crystal surface. This creates a thin (0.1- 
0.2 microns) alpha emitting layer. Although a high concentration of radioactive nuclei 
(0.1-1 ppm) is needed to account for the observed rate, the predicted spectrum agrees 
quite well with observations. Note that it is not known how such a large concentration 
of radioactive nuclei can appear on the surface of an encapsulated crystal. 

Similar hypotheses on the source of the anomalous fast events have been suggested 



by Saclay groups [O, 15] and by Cooper et al. |lb[. However, note that the hypothesis 



that anomalous events are due to 214 Po decay [16| requires a constant supply of radon 
because of the short lives of isotopes decaying into 214 Po. 

If the source of fast events is indeed on the surface of the crystal, then it can 
be removed by polishing the surface. This is hard to do with encapsulated Nal(Tl) 
crystals but such an experiment can be done with CsI(Tl) providing it shows a similar 
rate of fast events. The advantages of CsI(Tl) crystals are: a) they are only slightly 
hygroscopic and can be easily handled; b) they show better discrimination capability 
between electron and nuclear recoils || [17|- (Note that Csl is generally ruled out for 
dark matter searches due to high intrinsic background). 

3. Test with CsI(Tl) crystal 

Tests were performed with an 0.8 kg CsI(Tl) Harshaw crystal previously studied in 
the laboratory to evaluate its characteristics such as quenching factor for recoils and 
discrimination power, relevant to dark matter searches. The results have been reported 
elsewhere 1131. 
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The crystal was subsequently moved to the underground laboratory at Boulby and 
tested for background rate and anomalous fast events. In all tests we applied the stan- 
dard procedure of pulse shape analysis adopted by the UK Dark Matter Collaboration 
for Nal(Tl) dark matter detectors |2|, |6], [5], ^ . Pulses from both PMTs were integrated 
using a buffer circuit and then digitised using a LeCroy 9350A oscilloscope driven by 
a Macintosh computer running Labview-based data acquisition software. The digi- 
tised pulse shapes (10 /is digitisation time) were passed to the computer and stored on 
disk. Final analysis was performed on the sum of the pulses from the two PMTs. Our 
standard procedure of data analysis involves the fitting of a single exponential to each 
integrated pulse to obtain the index of the exponent, r. Although scintillation pulses 
from CsI(Tl) have an additional second component ||, the pulses can nevertheless be 
well fitted by a single exponential if fits are restricted to data below 1500 ns. This 
fraction of the pulse contains the major contribution to the integrated pulse amplitude 
so that the distortion of the fit due to the presence of the second exponential at large 
time scales was found to be insignificant ||. This approximation has the advantage 
that a three parameter fit can be used on each pulse and a simple discrimination pa- 
rameter defined, rather than a considerably more complicated six parameter fit in the 
case of two decay constants. The free fit parameters used are: the time constant of the 
single exponent, r; a normalisation constant and the start time of the pulse. 

For each run the distribution of the number of events versus the time constant of 
the exponent (r) was generated for a range of energy bins, r-distributions for each 
population of pulses can be approximated by a gaussian in ln(r) (log(Gauss) function) 
H, [| (for a more detailed discussion of the distributions see [[l(J and references therein): 
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The CsI(Tl) r-distributions were fitted with this gaussian in ln(r) with the three free 
parameters: time constant r , width w and normalisation factor A^ . In experiments 
where a second population is seen (for example, nuclear recoils from a neutron source or 
anomalous fast events), the resulting r-distribution can be fitted with two log(Gauss) 
functions with the same width w. 

In the low background conditions of the underground laboratory at Boulby, the 
CsI(Tl) crystal was found to show an anomalous population of fast events. Figure [5]a 
shows time constant distribution of events with visible energy 30-50 keV together with 
a fit to a log(Gauss) function (14.3 kg x days of exposure). The spectrum of these 
events presented in Figure ^ (crosses) has been traced up to MeV energies. The rate 
and shape of the spectrum are similar to those observed in the Nal(Tl) encapsulated 
detectors |T||. A peak at about 2.9 MeV corresponds to the 5.3 MeV alphas from 
210 Po — > 206 Pb a-decay, assuming a quenching factor for alphas of 0.55. (Scintillation 
efficiency of 0.50 ± 0.05 for alphas was measured for this crystal with 241 Am and 137 Cs 
sources). Note the absence of higher energy events which may be associated with the 
decay channels prior to 210 Po. This is in contradiction to what was suggested in Ref. 

After 2 months of running at various dynamic ranges the crystal was removed, 
polished and put into a sealed vessel with nitrogen atmosphere. After polishing the 
crystal was exposed to air for only a few hours during the installation procedure. 

The subsequent runs revealed a decrease in the rate of fast events by about a factor 
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of 4 (squares in Figure |]). The first two points below 100 keV show an upper limit 
to the rate. An accurate measurement of the rate at these energies is difficult because 
of the small mass of the crystal and the high rate of 7-background observed due to 
internal contamination of 137 Cs. The time constant distribution for events of 30-50 
keV after polishing is shown in Figure ^b (22.3 kg x days of exposure). It can be seen 
that the rate of anomalous fast events is significantly reduced (see Figure |3|a), though 
not completely suppressed probably due to the difficulty of removing the hard surface 
layer of Csl. 

17 high-energy events (visible energy of 5-6 MeV) were also detected during the first 
day after polishing. These are double-pulse events where the first pulse corresponds to 
the /3-decay of 212 Bi (half-life 1 hour) and the second due to the a-decay of 212 Po (half- 
life 0.3 fjis). These events are probably caused by contamination of the crystal surface 
during installation. No more of these events were seen after the first day (half-life of 
parent isotope 212 Pb is 10.6 hours). Visible energy of these double-pulse events agrees 
with the assumption that the scintillation efficiency of alphas is about 0.5. 

Only one prominent peak is seen in the spectra shown in Figure [|. The peak is 
probably due to 210 Po — > 206 Pb a-decay (5.3 MeV as). No decay chains 222 Rn — > 218 Po 
— ► 214 Pb or 224 Ra — > 220 Rn — > 216 Po — > 212 Pb were seen before or after polishing. This 
suggests that the concentration of U, Th and Ra in the bulk of the crystal is very low 
(less than 0.1 ppb). 

At least several months of exposure to Rn is needed to explain the rate of a-events 
in CsI(Tl) and Nal(Tl) detectors. This is not surprising for an unencapsulated Csl 
crystal but is hard to explain for Nal sealed detectors. 

4. Array of unencapsulated Nal(Tl) detectors - NAIAD 

The results obtained with the CsI(Tl) crystal at Boulby clearly indicate the im- 
portance of having access to the crystal surface for polishing and control. Such access 
can be granted by running unencapsulated crystals in high purity mineral oil or dry 
gas inside sealed plastic or copper vessels. Laboratory tests have shown also that high 
light yield (up to 10 photoelectrons per keV) can be reached with the aforementioned 
detectors [|, 0. 

Since 1999 the UKDMC has been developing a programme to run several unen- 
capsulated Nal(Tl) crystals mounted in an array - NAIAD (NAI Advanced Detector). 
The NAIAD array is designed to be flexible enough to allow various modes of operation 
with crystals. To date two types of module have been constructed: a "vertical" module 
filled with high purity mineral oil to protect crystals from moisture, and a "horizontal" 
module in which either an oil or dry nitrogen is used around the crystal. Operation of 
an encapsulated crystal is also possible in the horizontal module. Design details and 
predictions of sensitivity to WIMPs are given in Ref. M. 

The first vertical module of NAIAD has been running since February 2000. It 
contains a 14 cm diameter x 15 cm length crystal (termed DM74) with mass of about 
8.5 kg. The crystal was polished before installation. The total exposure (excluding 
calibration runs) is 996 kg x days. Figure [| shows a typical distribution of time 
constants for events with 35-40 keV visible energy together with a fit to a log(Gauss) 
function. As with the polished CsI(Tl) crystal the rate of anomalous fast events is 
greatly suppressed. The calculated limit on the rate of anomalous fast events as a 
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function of visible energy is shown in Figure ^| together with the energy spectra of fast 
events measured in typical encapsulated crystals. 

A second (horizontal) module containing 4 kg unencapsulated and polished crystal 
in nitrogen (DM72) also does not show the presence of fast events. This second detector 
is currently running underground at Boulby. Analysis of data from both crystals in 
terms of limits on the WIMP-nucleon and WIMP-proton cross-sections is in progress. 

5. Pulse shape analysis versus annual modulation 

Pulse shape analysis (PSA) is not the only technique used with Nal(Tl) detectors 
for dark matter searches. The DAMA group J12] searches for an annual modulation in 



the background counting rate of their Nal(Tl) array without PSA (PSA is used only 
to discriminate between scintillation pulses and PMT noise). Evidence for such an 



annual modulation in the background rate has been reported by DAMA (see fll|] and 
references therein) indicating a possible signal from WIMPs. 

It would appear that the DAMA group could, in fact, confirm or exclude this pos- 
sibility using pulse shape analysis. In Ref. DAMA presented results of the annual 



modulation analysis (positive signal) together with previous limits on WIMP-nucleon 
cross-section obtained with PSA. All sets of data (DAMA/NaI-0 analysed using PSA 
and DAMA/NaI-1 - DAMA/NaI-4 analysed using annual modulation) were obtained 
with the same experimental set-up and under similar conditions such as background 
rate etc. Pulse shape analysis applied to the first data set (DAMA/NaI-0) allowed 
DAMA to put limits on WIMP-nucleon interactions of the order of (5 — 6) x 
pb for 50-100 GeV mass WIMPs with halo density of 0.3 GeV/cm 3 (see JT2|, |13| for 
details and halo parameters used). This limit was obtained with 4123.2 kg x days 
exposure. The subsequent data sets (DAMA/NaI-1 - DAMA/NaI-4) totaling 57986 kg 
x days showed a positive signal at 4 a confidence level using annual modulation anal- 
ysis without PSA. WIMP parameters derived from this analysis are: My/ = (52+g°) 
GeV and cross-section a p = (7. 2^9) x 10~ 6 pb with the same halo parameters fl2|j . 
Simple statistical considerations show that PSA applied to this 15 times larger expo- 
sure from all five data sets (DAMA/NaI-0 - DAMA/NaI-4) compared to the first one 
(DAMA/NaI-0) could yield a limit on the WIMP-nucleon cross-section improved by a 
factor of 3.9. Such an analysis would allow DAMA either to confirm the modulated 
signal or to exclude practically the whole region of parameters that they derive from 



the observed modulated signal (see figures 3 and 4 in Ref. (T^j for the DAMA allowed 
region of WIMP parameters). The same statistical considerations suggest that the 
signal reported in Ref. [19| for the period DAMA/NaI-1 with 4549 kg x days should 



have been at a 1.1 a confidence level, to be consistent with the signal observed at 4 a 
confidence level with 57986 kg x days exposure. 



6. Conclusions 

Tests with a CsI(Tl) and unencapsulated Nal(Tl) crystals have shown that anoma- 
lous fast events seen in several Nal(Tl) detectors at Boulby were probably due to surface 
as. Radioactive a-emitting isotopes had likely been implanted into crystal surfaces by 
radon decay. Polishing the crystal surfaces removed a major part of the fast events. 
A new array of unencapsulated Nal(Tl) crystals (NAIAD) is being installed in the 
underground laboratory at Boulby. Data from the first of these modules do not reveal 
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anomalous fast events. 
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Figure 1: a) Time constant distribution for events with visible energy 35-40 keV from 
one encapsulated Nal(Tl) detector; b) similar distribution for Compton events from a 
gamma source. Solid curves show fits to Gaussian distributions on a logarithmic scale 
(log(Gauss)-function) . 
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Figure 2: a) Time constant distribution for events with visible energy 40-45 keV from 
encapsulated Nal(Tl) detector; b) similar distribution for calibration run with neutron 
source. Solid curves show fits to a sum of 2 log(Gauss)-functions. Dashed curves show 
fits to gamma-induced events. Fits to anomalous fast events (a) and neutron-induced 
events (b) are plotted by dotted curves. 
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Figure 3: a) Time constant distribution for events with visible energy 30-50 keV from 
the CsI(Tl) crystal; b) similar distribution after crystal polishing. Solid curves show 
fits to a log(Gauss)-function. 
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Figure 4: Rate of fast events (as) in CsI(Tl) crystal before (crosses) and after (squares) 
polishing. The first two points after polishing show limits at 90% confidence level. 
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Figure 5: Time constant distribution for events with visible energy 35-40 keV from the 
first NAIAD module (DM74). Solid curve shows a fit to a log(Gauss)-function. 
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Figure 6: Rate of anomalous fast events in Nal(Tl) encapsulated crystals and limits on 
the rate of these events in the unencapsulated crystal (DM74): triangles - UKDMC 
DM46 crystal; diamonds - Saclay crystal, also named DM70 (see also [11. 14]); arrows 
- limits at 90% confidence level from DM74 data. 
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